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ABSTRACT 

The BL Lac population is set apart from the rest of the active galactic nuclei by a number of peculiar 
features, including remarkably little signs of cosmological evolution. We focus on this feature, and 
take hint from the observations to reconsider extraction of the primary energy from a Kerr hole via 
the Blandford-Znajek (BZ) process; in fact, we stress how this can power the emitting jets of the BL 
Lac objects in a regime of low accretion. We show that the timescale for the BZ process to exhaust 
the rotational energy of the hole spans many gigayears in all such sources. We compute scales long 
and uniform enough to yield an intrinsically slow population evolution, giving rise to number counts of 
"Euclidean" or flatter shape; even longer timescales apply to the additional, analogous process involving 
the inner accretion disk directly. Such slow evolution takes place in ambients and at epochs (z 1) 
where interactions of the host galaxies with companions become rare or weak. On the other hand, at 
higher z, we argue that frequent interactions of the hosts within groups enhanced, each for some 1CU 1 
Gyr, the accretion rates and the emissions from the disks, attaining power levels that are typical of the 
radio-loud quasars. But such interactions were also likely to disorder the large-scale magnetic field in the 
disk and to decrease the angular momentum of the hole, thus limiting the lifetime of all BZ emissions to 
under a few gigayears. Then the stronger evolution of the flat-spectrum radio-loud quasars was driven, 
similar to the radio-quiet ones, by the interaction rate decreasing rapidly with z. 

Subject headings: BL Lacertae objects: general - galaxies: active - galaxies: evolution - radiation 
mechanisms: nonthermal 



1. INTRODUCTION 

The BL Lacertae objects among the active galactic nu- 
clei (AGNs) exhibit a number of peculiarities (see Urry & 
Padovani 1995, Sambruna 1999, and references therein), 
namely: (a) lack of emission lines with equivalent width 
exceeding 5 A; (6) a pure nonthermal continuum; (c) flat- 
spectrum radio emission, and powerful 7-rays that extend 
up to TeV energies in some objects; (d) rapid variability; 
and (e) strong and highly variable optical polarization. To 
some extent, they share the three latter features with the 
flat-spectrum radio-loud quasars (FSRQs). 

The current understanding links the peculiarities c, d, 
and e to synchrotron and inverse Compton emissions from 
a relativistic jet with bulk Lorentz factors F ~ 3 — 10, 
viewed at a small angle of order T _1 (see Ghisellini et al. 
1993). On the other hand, the peculiarities a and b are 
widely related to a dearth of gas surrounding the activity 
center at both parsec and 10 /ipc scales, and/or to the lack 
of an adequate nonbeamed continuum. An additional pe- 
culiarity is constituted by the behavior of the population, 
which shows little or no signs of evolution (see Stickel et 
al. 1991; Wolter et al. 1991; Bade et al. 1998), at striking 
variance with the rest of the AGNs. 

To us, the peculiarities a and b suggest that current ac- 
cretion is low in BL Lac objects, providing little gas and 
a weak ionizing continuum. Taking up previous proposals 
(see Blandford 1993 and references therein), we consider 
the main power source for the BL Lac emissions to be 
provided by the extraction in electromagnetic form of ro- 
tational energy, stockpiled in the central black hole and 



associated with the inner accretion disk. We develop this 
view to explain the other BL Lac peculiarities, and specif- 
ically the weak evolution. 

2. EVOLUTIONARY BEHAVIOR 

Once Doppler enhancements of order r 3 ~ 4 are ac- 
counted for, the observations of BL Lac objects (see Sam- 
bruna, Maraschi & Urry 1996), indicate typical rest frame 
luminosities L $ 10 44 ergs s , which are moderate on the 
quasar standards; meanwhile, the rotational energies ex- 
tractable from Kerr holes (KHs) are around E ~ 5 10 61 M§ 
ergs. So their exhaustion would take in individual ob- 
jects times t ~ E/L of many gigayears, which are quite 
longer than the flashes lasting about 10 _1 Gyr of stan- 
dard quasars fueled by accretion (see Cavaliere & Vittorini 
1998). In the few instances of larger L, the Doppler factors 
inferred from the data in the framework of the synchrotron 
self-Compton radiation - which constitute lower bounds 
anyway - are particularly low and uncertain pending con- 
firmation from superluminal velocities (P. Padovani 1998, 
private communication). The long debated kinetic com- 
ponent of the jets is unlikely to increase the power budget 
of the typical BL Lac objects beyond some 10 45 ergs s~ 4 
(Ghisellini & Celotti 1998). 

The population-averaged (L(t)) is constituted by the 
convolution of the individual behaviors. When these are 
all similar, as we shall show below, (L(t)) also shares the 
same timescale r. Heuristically, a long r will imply little 
population evolution and, therefore, relatively flat number 



counts, or equivalently V/V ma x close to 0.5. 
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A quantitative statement can be made in simple terms 
on recalling that the integral counts N(> S) at high fluxes 
S may be cast into the form (see Weinberg 1972): 



L h oc B 



A 2 for the power extracted from the KH it- 



self in the form of Poynting flux; here A < 1 denotes the 
angular momentum of the hole, J, in units of GM 2 /c, 
and B is the magnetic field threading the horizon at 



N(> S)ocS-V 2 [l-C(S o /S) 1 / 2 + 0(S- 1 )] . (1) r h = [1 + (I - A 2 ))^ 2 ]GM/c 2 



In turn, C can be expressed (Cavaliere & Maccacaro 1990) 
in terms of the time scale r e for population changes by 



C = 



3D <1 2 > 
4R H < P/ 2 > 



[2(1 + a) . (2) 

tloTe 



Here Rh — 3000 h 1 Mpc is the Hubble radius with 
the Hubble constant h = ff o /100 km s" 1 Mpc -1 , D Q = 
{P /AttS ) 1/2 « 200 h' 1 Mpc in terms of the power and 
flux of the brightest BL Lac objects, < l n > is the (normal- 
ized) n-th moment of the luminosity function (LF), and a 
is the spectral index. 
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FIG. 1 - The runs of A(t) and of L(t) for M = const, correspond- 
ing to initial normalized accretion rates rn = 5 10~ 2 (upper curves), 
10~ 2 (middle curves), and 10 -3 (lower curves). Note the convergent 
behaviors at the lowest rates. 

The evolutionary rate r" 1 = t^, 1 + ([3 — 1)777 1 comprises 
two contributions that are due to density and to luminosity 
evolutions, respectively; the latter is amplified by a steep 
LF (modeled with L^ 13 in the relevant range), and td is 
negative in case of negative evolution. From cq. (2) it is 
seen that the bright counts tend to the "Euclidean" shape 
N(> S) cx S~ 3//2 or flatter even when (3 ^ 2.5, and even 
in the presence of positive evolution, provided this is slow 
enough to statisfy r e H Q > [2(1 + a)] -1 . 

3. A SIMPLE MODEL: BEAMS FROM KERR HOLES 

We begin with considering the specific model where the 
primary power is extracted from the rotational energy 
of a KH, following Blandford & Znajek (1977). They 
tackled a stationary, axisymmetric and force-free mag- 
netosphere surrounding the KH, and derived the scaling 



The strength of the field that is sustainable at the hori- 
zon constitutes a delicate issue. A recent reappraisal by 
Ghosh & Abramowicz (1997) has been based on continu- 
ity of the magnetic pressure at the inner edge of the ac- 
cretion disk and on numerical computations - which are 
self-consistent, if with dynamic range limited to small and 
medium scales - of its MHD and thermal structure. They 
find that B 2 /8it ~ ap holds for standard a-disks in terms 
of the maximal pressure at the inner rim; they (conserva- 
tively) evaluate the resulting power to read 



2 1 44 M 8 yl 2 ergs s" 1 

8 lO^Mg 11710 ™^ 2 ergs S" 1 , 



(3) 



when the inner region is dominated, respectively, by 
radiation [p cx (aM)^ 1 ] or by gas pressure [p oc 
(aM)- 9 / 10 m 4 / 5 ]. We recall from Novikov & Thorne 
(1973) that only when the accretion rate in Eddington 

units exceeds m = Mc 2 /L E ~ 310- 4 M 8 ~ 1/8 (for A = 1) 
does the transitional radius oc M 16 / 21 (separating the 
radiation-dominated region from the gas-dominated re- 
gion) exceed the radius r ms of the marginally stable orbit 
so that an inner, radiation pressure dominated region in 
fact occurs. Then the ratio of the Blandford- Znajek (BZ) 
emission to the thermal luminosity scales as Lh/L t h ~ 
310~ 2 A 2 /to. Thus, in the range 10~ 3 < m < 10~ 2 , the 
rest frame power Lh is adequate to sustain the nonther- 
mal, beamed emissions of BL Lac objects and to match or 
dominate the thermal component L t h- 

As for the evolutionary behavior of the population, we 
first consider the time dependence of Lh(t) for an individ- 
ual source; this is related, by equations (3), to the behav- 
iors of M(t) and of J(t). The key point is that accretion of 
J also takes place along with M , and is always important, 
as first pointed out by Bardeen (1970), being related to M 
by the conditions in the disk at r ms . 

In fact, the run of M(t) is described by 



dM 



e ms M - L h /c 2 



(4) 



where e ms ~ 1 is the adimensional specific energy at r ms ; 
actually, e ms = 0.58 holds for A = 1, rising to 0.94 for 
A = (dc Felice 1968; Bardeen 1970). On the other hand, 
the run of J(t) is coupled to the former by 



dJ ■ ,7 

— = JmsM 

at 



k Qt, 



(5) 



(Malquori 1996; Moderski & Sikora 1996). Here j ms is 
the specific angular momentum at r ms , with values close 
to 2 e m sCrh\ moreover, Vlh = Ac/2rh, while k — 1/2 for 
maximum power extraction. 

On the right-hand side of equation (5) the accretion of 
J can be balanced out by the loss due to the BZ pro- 
cess, causing dJ/dt to vanish for values of A around 50 to, 
given 10" 3 < to < 10~ 2 . Instead, A(t) oc J{t)/M 2 {t) de- 
creases monotonically (see Fig. 1). Correspondingly, Lh{t) 
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(shown in Fig. 1 with all relativistic factors and detailed 
yl-dcpcndcncies in place) decreases steeply at first because 
of the relativistic factors, but then it tends to level off (or 
to increase again for high M) because of the hole mass 
increase; in between, the decrease takes place on the scale 
r = L/L ~ J/2 J. Such behavior and scale are similar for 
all objects at low m where the BZ power dominates. 

In sum, the sources powered by KHs fade slowly, but re- 
main above levels L/jo/10 ~ 10 44 erg s _1 that are relevant 
for BL Lac objects for times r ^ 8 Gyr. So their popula- 
tion undergoes a slow luminosity evolution; the condition 
r e H a ~ [2(1 + a)] -1 tends to hold especially for h ~ 0.7, 
and the bright counts in any band ought to be quite flatter 
than with the quasars. 

As an example, we compute in detail eqs. (1) and (2) 
for the radio band on the basis of the above luminosity 
evolution alone. In the radio, the BL Lac counts are ar- 
guably least sensitive to the spectral changes with power 
proposed by Fossati et al. (1998; see also Georgenopoulos 
& Marscher 1998 and Rector et al. 1999). There a ~ 
and (3 ~ 2.5 apply (we refer to the model LF computed by 
Padovani & Giommi 1995), and the condition r e H ~ 1/2 
for nearly Euclidean counts is closely met. In fact, we find 
the shape shown in Figure 2 by the solid line; this is easily 
consistent with the best fit {dashed line) to the counts by 
Stickel et. al. (1991), a published sample that is com- 
plete although it is affected by much statistical variance. 
A flatter LF will yield counts flatter still. On the other 
hand, a flatter LF, but with fast evolution, would yield 
steep counts as illustrated by the dotted line in Figure 2. 

Actually, our results provide an upper limit to the 
counts, since the values of r evaluated above provide min- 
imal timescales while an additional negative density evo- 
lution is likely to occur, as discussed next. More extensive 
and deeper data will be released soon (see Giommi, Menna 
& Padovani 1999 and Rector et al. 1999) and will provide 
more challenging evolution indicators for our prediction. 

4. DISCUSSION AND CONCLUSIONS 

All peculiarities of the BL Lac objects are consistent 
with the extraction in electromagnetic form of the rota- 
tional energy stockpiled in a KH, in a regime of low ac- 
cretion. In fact, a primary Poynting flux, as produced via 
the BZ process by the KH itself, is attractive (Ghiscllini 
& Celotti 1998) because it is energetically adequate and is 
conducive to jets with bulk L ~ 10; by way of the implied 
large-scale electric fields, it is conducive also to distributed 
and continuous acceleration up to isotropic electron ener- 
gies in the range of 10 2 GeV. Thus, highly beamed, non- 
thermal (and polarized) continua extending into the TeV 
7-ray range can be generated; if a jet is oriented close to 
our line of sight, high apparent power and fast time vari- 
ability will be observed (Bcgclman et al. 1984). In turn, 
low current accretion rates, m ^5 10 -2 , are consistent with 
negligible dilution by thermal continua (then also unable 
to degrade the TeV 7-rays) and with emission lines of small 
equivalent width. But our thrust has been to show how 
this view leads us to predict intrinsically slow population 
evolution. 

Actually, an additional MHD emission Ld is produced 
by the inner disk (Blandford & Znajek 1977; Livio, Ogilvie 
& Pringle 1999). This scales still with B 2 , but also with 
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rl^rj/ 2 so that the ratio Ld/Lh oc (r^/r^) 3 / 2 cx M 1 - 15 ; the 
coefficient is independent of A 2 (t) and somewhat larger 
than in eq. (3a). Consideration of Ld does not affect 
materially our evaluations as long as m 10~ 2 , and so 
Ld ~ Lh, which is appropriate to the BL Lac objects. 
Then the KH and the disk are dynamically and magnet- 
ically linked nearly in one system in which - within (un- 
certain) factors of a few - the hole still may be regarded 
as the energetic flywheel. 
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FIG. 2 - BL Lac counts at 5 GHz, from Stickel et al. 1991. The 
dashed line reproduces their best fit N(> S) oc S -1,6 (the actual 
slope is 1.58 ± 0.44 at a Iff level). The solid line gives our result 
for pure luminosity evolution with r = 8 Gyr, assuming constant 
Lradiol Lh and P = 2.5. The dotted line represents for comparison 
a case where (i = 2.2 and r = 2 Gyr. 

The stand in this context of the other Blazars, the FS- 
RQs, is interesting as such objects exhibit features similar 
to the BL Lacs, less pure but with larger outputs. Tap- 
ping the rotational energy in order to produce relativistic 
jets, blazing beams, and nonthermal luminosities is tanta- 
lizing also for these objects, but there are also important 
differences. 

Now the maximal KH power Lh from eq. (3a) is inad- 
equate - with masses limited to M 410 9 h~ x M & (see 
Magorrian et al. 1998) - to match powers of 10 46 h~ 2 erg 
s~ 4 or more; a larger disk component Ld to the BZ emis- 
sion is cogently called for. In fact, Ld 3> Lh now holds 
(with the KH now slaved to the disk) since m <~ 1; we 
know this is to hold for the FSRQs, because their emission 
lines and thermal components, which dilute the BL Lac- 
like peculiarities to approach the standard quasar activity, 
do signal currently high accretion rates. We also expect 
the larger disk component in FSRQs to imply - relative to 
the BL Lac conditions - more plasma lifted up and elec- 
tric fields more accurately screened out; then the electrons 
not only may loose their energy faster at such larger total 
outputs (Fossati et al. 1998), but are also accelerated to 
lower energies, which is another, intrinsic reason to emit 
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softer 7-rays. 

Since the BZ emissions are marked by long timescales at 
constant M, the long persistence of radio-emitting, rota- 
tionally driven jets might be seen as statistically embarass- 
ing (see Moderski, Sikora & Lasota 1998). However, we 
argue that the actual timescales of the FSRQs are much 
shorter because of the boundary conditions set to such 
high values of M by the host galaxy For the sake of this 
argument, we shall discuss first the scenario (most con- 
straining to us) where the BL Lac objects are descendants 
of FSRQs. 

We recall from Cavaliere & Vittorini (1998) that recur- 
rent interactions and merging events of the E host galaxies 
with satellites or with group companions (for which there 
is plenty of morphological and kinematical evidence; see 
Bertola & Corsini 1998) are necessary to cause substan- 
tial losses of angular momentum, |Aj'|/j ~ 1, of the gas 
residing in either interaction partner; such events every 
1 — 2 Gyr trigger for some 10 _1 Gyr gas inflow toward the 
nucleus, which is conducive to M and also to accretion of 
angular momentum onto the hole. We concur with Moder- 
ski, Sikora & Lasota (1998) that the sign of such imported 
momentum relative to the previous one makes a difference 
as for disk structure and for the KH rotation. We stress 
that each strong interaction has a 1/2 chance to accrete 
a component opposite to the initial J of the KH and of 
the disk, thus wrecking the latter's large-scale magnetic 
field and spinning the former down, effectively quenching 
their emissions. When the disk is rebuilt, the internal con- 
ditions for large-scale B and those in the host galaxy for 
beam propagation (see Livio 1999) may not be met again. 

Thus we hold that the probability is small for combining 
persistent extraction of BZ power with frequent episodes of 
high m. But the quasar evolution since z *=s 3 is governed 
just by high-m episodes driven by interactions; these are 
initially frequent and strong in the era of group formation, 
but then drop off in rate and strength with decreasing z as 
groups are reshuffled into clusters (Cavaliere & Vittorini 
1998). As long as such interactions arc effective, they will 
intermittently switch off all BZ outputs from the KH and 
the disk. 

We conclude that strong evolution and long-lived extrac- 
tion of BZ power entertain a complementarity relationship. 
Indeed, the rotational power and the beams it feeds persist 
safely for many gigayears only in a regime of permanently 
low accretion, which is what prevails for z 1.5, after a 
"last interaction" setting A ~ 1; such conditions will lead 



to the BL Lac phenomenon. At higher z, the most likely 
activity pattern was that corresponding to FSRQs, with 
dominant thermal and disk BZ emissions limited to a few 
gigayears by initially frequent, but then subsiding, inter- 
actions; the evolution of such sources is then expected to 
be close (cf. Goldschmidt et al. 1999) to the standard 
quasars', being driven by the same z-dependent (m) only 
somewhat slowed down by the beam attempts to persist 
in between interactions. On the other hand, if it is main- 
tained (see Jackson & Wall 1999) that BL Lacs and FSRQs 
belong to different and secluded lineages, then little con- 
straints will arise from FSRQs for the the present view of 
the astrophysics and evolution of BL Lac objects. 

Finally, enlarging our view of this recently rejuvenated 
field above details and beyond the theoretical dream of 
a relativistic homopolar inductor, we summarize and em- 
phasize the two key issues in this context. First and con- 
cerning the BL Lac population, up to now flat counts and 
V/V max ~ 0.5 persistently signal a weaker evolution than 
with the quasars. This indicates either a slower luminosity 
evolution with timescales tl/({3 — 1) of several gigayears, 
as considered above, or a negative density evolution with 
to < 0, or some mix of the two (some negative density 
evolution is also implied from BL Lac birth since z « 1.5). 
Second, concerning the source structure, the current wis- 
dom on disks (but ADAF or ADIOS flows, possibly set- 
ting in at low rh, as discussed by Blandford & Begelman 
1998, have yet to be firmly fit in) has that the disk ra- 
diation pressure is independent of accretion for m above 
some 10 ~ 3 , since then p oc T 4 holds in the inner region. 
As long as good coupling of radiation, particles and large- 
scale magnetic fields can be retained, this leads to a large 
ratio (Lh + Ld)/L t h oc 1/m; in addition, as m approches 
10 -3 from above, shrinks to approach r ms , and com- 
parable L h ~ Ld obtain. Such conditions are best suited 
to yield pure nonthermal radiation of power moderate but 
extending up to TeV 7-rays, and by the same token they 
yield persistent beams with intrinsically weak evolution, 
precisely the features marking out the BL Lac objects. 
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